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We present the first relativistic calculations of the final phase of inspiral of a binary system consist- 
ing of two stars built predominantely of strange quark matter (strange quark stars). We study the 
precoalescing stage within the Isenberg- Wilson-Mathews approximation of general relativity using 
a multidomain spectral method. A hydrodynamical treatment is performed under the assumption 
that the flow is either rigidly rotating or irrotational, taking into account the finite density at the 
stellar surface — a distinctive feature with respect to the neutron star case. The gravitational- 
radiation driven evolution of the binary system is approximated by a sequence of quasi-equilibrium 
configurations at fixed baryon number and decreasing separation. We find that the innermost sta- 
ble circular orbit (ISCO) is given by an orbital instability both for synchronized and irrotational 
systems. This constrasts with neutron stars for which the ISCO is given by the mass-shedding limit 
in the irrotational case. The gravitational wave frequency at the ISCO, which marks the end of 
the inspiral phase, is found to be ~ 1400 Hz for two irrotational 1.35 Mq strange stars and for the 
MIT bag model of strange matter with massless quarks and a bag constant B = 60 MeVfrn -3 . 
Detailed comparisons with binary neutrons star models, as well as with third order Post-Newtonian 
point-mass binaries are given. 

PACS numbers: 04.40.Dg, 04.30.Db, 04.25.Dm, 97.10.Kc, 97.60. Jd 



I. INTRODUCTION 

One of the most important prediction of general rel- 
ativity is gravitational radiation. Coalescing neutron 
star binaries are considered among the strongest and 
most likely sources of gravitational waves to be seen by 
VIRGO /LIGO interferometers |2( . Due to the emis- 
sion of gravitational radiation, binary neutron stars de- 
crease their orbital separation and finally merge. Gravi- 
tational waves emitted during the last few orbits of inspi- 
ral could yield important informations about the equa- 
tion of state (EOS) of dense matter [1 H, H . With 
accurate templates of gravitational waves from coalesc- 
ing binary compact stars, it may be possible to extract 
information about physics of neutron stars from signals 
observed by the interferometers and to solve one of the 
central but also most complex problem of physics — the 
problem of the absolute ground state of matter at high 
densities. It is still an open question whether the core of 
a neutron star consists mainly of superfluid neutrons or 
exotic matter like strange quark matter, pions or kaons 
condensates (see e.g. Ref. @| for a recent review). The 
possibility of the existence of quark matter dates back 
to the early seventies. Bodmer || remarked that mat- 
ter consisting of deconfined up, down and strange quarks 
could be the absolute ground state of matter at zero pres- 
sure and temperature. If this is true then objects made 
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of such matter, the so-called strange stars, could exist 
0, El [HI ■ Strange quark stars are currently considered 
as a possible alternative to neutron stars as compact ob- 
jects (see e.g. 0,0,0 and references therein). 

The evolution of a binary system of compact objects 
is entirely driven by gravitational radiation and can be 
roughly divided into three phases : point-like inspiral, hy- 
drodynamical inspiral and merger. The first phase corre- 
sponds to large orbital separation (much larger than the 
neutron star radius) and can be treated analytically using 
the post-Newtonian (PN) approximation to general rela- 
tivity (see Ref. 0] for a review). In the second phase the 
orbital separation becomes only a few times larger than 
the radius of the star, so the effects of tidal deformation, 
finite size and hydrodynamics play an important role. In 
this phase, since the shrinking time of the orbital radius 
due to the emission of gravitational waves is still larger 
than the orbital period, it is possible to approximate the 
state as quasi-equilibrium [0| ■ The final phase of the 
evolution is the me rger of the two objects, which occur 
dynamically [l8i flflt I2CM . f2ll | . Note that quasi-equilibrium 
computations from the second phase provide valuable ini- 
tial data for the merger H GJ, IU E| 

Almost all studies of the final phase of the inspiral 
of close binary neutron star systems employ a simpli- 
fied EOS of dense matter, namely a polytropic EOS 
0, II E3, El ES ISII2E ill23, El There are only two 
exceptions: (i) Oechslin et al. have used a pure nuclear 
matter EOS, based on a relativistic mean field model and 
a 'hybrid' EOS with a phase transition to quark mat- 
ter at high density @; (ii) Bejger et al. have computed 
quasi-equilibrium sequences based on three nuclear mat- 
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ter EOS In this article we present results on the hy- 
drodynamical phase of inspiraling binary strange quark 
stars described by MIT bag model. The calculations are 
performed in the framework of Isenberg- Wilson- Mathews 
approximation to general relativity (see Ref. |29( for a 
review). We consider binary systems consisting of two 
identical stars. We choose the gravitational mass of each 
star to be 1.35 Mq in infinite separation in order to be 
consistent with recent population synthesis calculations 
[30| and with the current set of well-measured neutron 
star masses in relativistic binary radio pulsars 012 . 
We compare the evolution of a strange star binary system 
with a neutron star binary in order to find any character- 
istic features in the gravitational waveform that will help 
to distinguish between strange stars and neutron stars. 
We consider two limiting cases of velocity flow in stellar 
interior: the irrotational and the synchronized case in or- 
der to exhibit the differences between these two extreme 
states. The irrotational case is more realistic since the 
viscosity of neutron star matter (or strange star matter) 
is far too low to ensure synchronization during the late 
stage of the inspiral |33l l34| . Due to the finite density 
at the surface of bare strange stars, we had to introduce 
a treatment of the boundary condition for the velocity 
potential (in the irrotational case) different from that of 
neutron stars, where the density vanishes at the stellar 
surface. 

The paper is organized in the following way: Sec. II 
is a brief summary of the assumptions upon which this 
work is based, Sec. Ill is devoted to the description of the 
EOS used to describe strange stars and neutron stars. In 
Sec. IV we briefly describe the basic equations for quasi- 
equilibrium and derive the boundary condition required 
for solving the fluid equation of irrotational flow with fi- 
nite surface density, which is relevant for strange stars. 
In Sec. V we present the numerical results for corotating 
and irrotational strange stars binaries and compare their 
quasistationary evolution with that of neutron stars, as 
well as with that of post-Newtonian point-masses. Sec- 
tion VI contains the final discussion. Throughout the 
paper, we use geometrized units, for which G = c = 1, 
where G and c denote the gravitational constant and 
speed of light respectively. 



II. ASSUMPTIONS 

The first assumption regards the matter stress-energy 
tensor T, which we assume to have the perfect fluid 
form: 



(e + p)u< 



pg, 



(i) 



where e, p, u and g are respectively the fluid proper 
energy density, the fluid pressure, the fluid 4- velocity, and 
the spacetime metric. This is a very good approximation 
for neutron star matter or strange star matter. 

The last orbits of inspiraling binary compact stars 
can be studied in the quasi-equilibrium approximation. 



Under this assumption the evolution of a system is ap- 
proximated by a sequence of exactly circular orbits. This 
assumption results from the fact that the time evolution 
of an orbit is still much larger than the orbital period 
and that the gravitational radiation circularizes an orbit 
of a binary system. This implies a continuous spacetime 
symmetry, called helical symmetry [35l l3r| represented 
by the Killing vector: 



dt dip' 



(2) 



where is the orbital angular velocity and d/dt and 
d/dip are the natural frame vectors associated with the 
time coordinate t and the azimuthal coordinate ip of an 
asymptotic inertial observer. 

One can then introduce the shift vector B of co- 
orbiting coordinates by means of the orthogonal decom- 
position of £ with respect to the E t foliation of the stan- 
dard 3+1 formalism: 



iVn-B, 



(3) 



where n is the unit future directed vector normal to St, 
iV is called the lapse function and n ■ B = 0. 

We also assume that the spatial part of the metric 
(i.e. the metric induced by g on each hypersurface S t ) 
is conformally flat, which corresponds to the Isenberg- 
Wilson-Mathews (IWM) approximation to general rel- 
ativity OH, |S3 (see Ref. [36] for a discussion). Thanks 
to this approximation we have to solve only five of the 
ten Einstein equations. In the IWM approximation, the 
spacetime metric takes the form: 



ds 2 



-(N 2 - BiB l )dt 2 - 2B t dt dx l + A 2 f tj dx l dx j , (4) 



where A is some conformal factor, /y the flat spatial 
metric and Latin indices run in {1, 2, 3} (spatial indices). 
The comparison between the IWM results presented here 
and the non-conformally flat ones will be performed in a 
future article |4*cj . 

The fourth assumption concerns the fluid motion in- 
side each star. We consider two limiting cases: synchro- 
nized (also called corotating) motion and irrotational 
flow (assuming that the fluid has zero vorticity in the 
inertial frame). The latter state is more realistic. 

We consider only equal- mass binaries consisting of 
identical stars with gravitational masses M\ = Mi = 
1.35 Mq measured in infinite separation. The main rea- 
son for choosing these particular masses is that five out of 
six observed binary radio pulsars have mass ratio close to 
unit y an d gravitational masses of each star ~ 1.3— 1.4M Q 
[31J, l32| . In addition population synthesis calculations 
30] have shown that a significant fraction of the observed 
binary neutron stars in gravitational waves will contain 
stars with equal masses ~ 1.4 Mq and systems consisting 
of a low and a high mass neutron star. 
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FIG. 1: Gravitational mass M versus areal radius R for se- 
quences of static strange quark stars described by the sim- 
plest MIT bag model (solid line) and neutron stars described 
by polytropic EOS with 7 = 2.5 and k — 0.0093 mo n^c 
(dashed line). The two sequences are crossing at the point 
M = 1.35 M Q and R = 10.677 km (marked by a circle). 



III. THE EQUATION OF STATE AND 
STELLAR MODELS 

It has been shown 0, Q that the evolution of equal- 
mass binary neutron stars depend mainly on the com- 
pactness parameter M/R, where M and R are the grav- 
itational mass measured by an observer at infinity for a 
single isolated neutron star and the stellar radius respec- 
tively. It is therefore interesting to check if the properties 
of inspiraling strange stars can be predicted by studying 
binaries consisting of polytropic neutron stars having the 
same mass and the same compactness parameter. There- 
fore we perform calculations for two different equations 
of state of dense matter: a strange quark matter EOS 
and a polytropic EOS. 

Typically, strange stars are modeled 0, 0] with an 
equation of state based on the MIT-bag model of quark 
matter, in which quark confinement is described by an 
energy term proportional to the volume |^3. The equa- 
tion of state is given by the simple formula 



p = a(p-p ), 



n(p) = n 



a p 



Po 



1/(1+0) 



(5) 



(6) 



where n is the baryon density and a, po,n>o are some 
constants depending on the 3 parameters of the model 
(the bag constant B, the mass of the strange quarks m s 
and the strenght of the QCD coupling constant a). In 
general this equation corresponds to a self-bound matter 
with mass density po and baryon density no at zero pres- 
sure and with a fixed sound velocity {\fa) at all pressures. 



It was shown that different strange quark models can be 
approximated very well by Eqs © and © 0> 113 ■ 

In the numerical calculations reported in the present 
paper we describe strange quark matter using the sim- 
plest MIT bag model (with massless and non-interacting 
quarks), for which the formula (J5J is exact. We choose the 
value of the bag constant to be B = 60 MeV fm -3 . For 
this model we have a = 1/3, po = 4.2785 x 10 14 g cm -3 
and no — 0.28665 fm~ 3 . In general for the MIT bag 
model the density of strange quark matter at zero pres- 
sure is in the range ~ 3 x 10 14 — 6.5 x 10 14 g cm~ 3 and 
a between 0.289 and 1/3 (for < a < 0.6 and < m s < 
250 MeV) ji^l ■ The higher value of a and of p the higher 
compactness parameter of a star with fixed gravitational 
mass. 

Up to now, the majority of calculations of the hydrody- 
namical inspiral phase H H liH I2I I21 153 EEi IM Wi pi 

and all calculations of the merger phase [lg, EH |2(J, ULLi 

have been performed for binary systems containing neu- 
tron stars described by a polytropic EOS: 



P 



(7) 



where k and 7 coefficients are some constant numbers: 
k represents the overall compressibility of matter while 
7 measures the stiffness of the EOS. The total energy 
density is related to the baryon density by 



e(n) 



7- 



1 



u n 



(8) 



where fio is the chemical potential at zero pressure. 

In order to compare results for strange stars with those 
for neutron stars, we determine the values of k and 7 
which yield to the same radius for the gravitational mass 
AI = 1.35 Mq as that obtained for a static strange star. 
It was shown |6j that the properties of inspiraling neu- 
tron stars described by realistic EOS can be, in a good 
approximation, predicted by studying binaries with as- 
sumed polytropic EOSs with 7 = 2 or 2.5. For a 1.35 
Mq strange star we have a high value of compactness 
parameter M/R = 0.1867 so we have choosen 7 = 2.5, 
for which we found k = 0.00937 monj^, with the rest 
mass of relativistic particles mo := 1.66 x 10~ 27 kg and 
"nuc = 0.1 fm~ 3 . 

In Fig. ^ we present the mass-radius relation for a se- 
quence of static stars described by the simplest MIT bag 
model (solid line) and the polytropic EOS (dashed line) 
parametrized by central density. For small mass strange 
stars M ~ R 3 since density is almost constant inside 
a star ~ po- In the top panel of Fig. [21 we show the 
mass density distribution inside the strange star (solid 
line) and the neutron star described by polytropic EOS 
(dashed line) having gravitational mass 1.35 Mq and 
areal radius 10.667 km (the configurations correspond- 
ing to the crossing point on Fig. P). The huge density 
jump at the surface of the strange star corresponds to 
po = 4B. The value of density at the surface describes 
strongly or weakly bound strange matter, which in each 
case must be absolutely stable with respect to 56 Fe. 
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FIG. 2: Mass density (top panel) and the adiabatic index 
7 (bottom panel) versus the radial coordinate r for a static 
strange quark star (solid line) and a polytropic neutron star 
(dashed line), having both a gravitational mass of 1.35 Mq 
and an areal radius R — 10.677 km (resulting in the com- 
pactness parameter M/R = 0.1867). The vertical dotted line 
corresponds to the stellar surface. 



of strange stars are very large and for p — > po we have 
7 = a + p/ (p — po) —>■ oo. The EOS of neutron stars 
for densities lower than ~ 10 14 g cm -3 (the crust) is well 
established Q . In the outer crust of an ordinary neutron 
star the pressure is dominated by the ultra-relativistic 
electron gas, so we have 7 = 4/3. The values of the local 
adiabatic index in the inner crust of a neutron star de- 
pends strongly on density and varies from 7 ~ 0.5 near 
the neutron drip point to 7 ~ 1.6 in the bottom layers 
near the crust-core interface. 

In our calculations we use equation of state in the form: 

n = n(H), e = e(H), p = p(H), (10) 

where H is pseudo-enthalpy (the log-enthalpy) defined 
l>v: 



H(n) := In 



P 



nEn 



(11) 



where the energy per unit baryon number is Eq = m 
for a polytropic EOS, and E = p /n = 837.26 MeV for 
strange quark model described above. For our model of 
strange quark matter we have: 



p = p (3c 4W +l)/4, p = Po (e 4 "-l)/4, 



IV. EQUATIONS TO BE SOLVED 



n e 



(12) 



We refer the reader to Ref. [ig] for the derivation of 
the equations describing quasi-equilibrium binary stars 
within the IWM approximation to general relativity. Af- 
ter recalling these equations, we mainly concentrate on 
the equation for the velocity potential of irrotational 
flows. Actually this equation has a different structure for 
strange stars than for neutron stars. This results from 
the non- vanishing of the density at the stellar surfaces of 
strange stars (cf. the top panel in Fig. EJ). 



A. The gravitational field equations 



An important quantity relevant for evolution of binary 
compact stars is the adiabatic index: 

7 = dlnp/dlnn. (9) 

We assume that matter is catalized so the adiabatic in- 
dex can be calculated directly from EOS (see Refs. 0] 
and [45j for discussion on different kind of adiabatic in- 
dices and corresponding timescales). Note that for the 
polytropic EOS given by Eq. J7J) the index 7 coincides 
with the adiabatic index of a relativistic isentropic fluid. 
Dependence of the adiabatic index 7 on stellar radii for 
both EOS is shown in the bottom pannel of Fig. [2J The 
adiabatic index of strange matter is qualitatively dif- 
ferent from the adiabatic index for polytropic EOS or 
for realistic EOS. The values of 7 in the outer layers 



The gravitational field equations have been obtained 
within the 3+1 decomposition of the Einstein's equa- 
tions 153, j taking into account the helical symmetry 
of spacetime. The trace of the spatial part of the Einstein 
equations combined with the Hamiltonian constraint re- 
sults in two equations: 

Af = AnA 2 {E + S) + A 2 K l3 K l i - V>vV, (13) 

A/3 = 4ttA 2 S + -A 2 K il K ij 
4 

-^(V^VV + V^V*/?), (14) 

where Vj stands for the covariant derivative associated 
with the flat 3-metric and A := V*Vj for the asso- 
ciated Laplacian operator. The quantities v and (3 are 
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defined by v := \nN and (3 :— ln(AN), and denotes 
the extrinsic curvature tensor of the t — const hypersur- 
faces. E and S are respectively the matter energy density 
and the trace of the stress tensor, both as measured by 
the observer whose 4- velocity is (Eulerian observer): 

(15) 
(16) 



E := T^n", 



S := .i 2 /'-"/; ; . 

In addition, we have also to solve the momentum con- 
straint, which writes 

1 



A7V l + ^V (VjN 3 ) 



-~16irNA 2 (E +p)U i 
+2NA 2 K ij Vj{3f3 - 4z/),(17) 



where N l := B l + il(d/dip) 1 denotes the shift vector of 
nonrotating coordinates, and U l is the fluid 3-velocity. 

B. The fluid equations 

Apart from the gravitational held equations, we have 
to solve the fluid equations. The equations governing the 
quasi-equilibrium state are the relativistic Euler equation 
and the equation of baryon number conservation. Both 
cases of irrotational and synchronized motions admit a 
first integral of the relativistic Euler equation: 



H + v - In T + In T = const. 



(18) 



where Tq is the Lorentz factor between the co-orbiting 
observer and the Eulerian observer and T is the Lorentz 
factor between the fluid and the co-orbiting observers 
(r = 1 for synchronized binaries). 

For a synchronized motion, the equation of baryon 
number conservation is trivially satisfied, whereas for an 
irrotational flow, it is written as 

= A 2 hr n U^ViH + £H 

x[V i ^W i (H-P)+A 2 hUlW i T n ], (19) 

where 'J is the velocity potential, h := exp(iJ), £ the 
thermodynamical coefficient: 

dkiH 



d\nn 

and T n denotes the Lorentz factor between the fluid and 
the Eulerian observer and Ufa is the orbital 3-velocity with 
respect to the Eulerian observers: 



(21) 



The fluid 3-velocity U l with respect to the Eulerian 
observer is equal to Uq for synchronized binary systems, 
whereas 



IP 



1 



A 2 T„h 



V * 



(22) 



for irrotational ones. 



C. Boundary condition for the velocity potential 

The method of solving the elliptic equation l|19|) for 
the velocity potential is different for neutron stars and 
strange stars. In the case of neutron stars, the coefficient 
in front of the Laplacian vanishes at the surface of 
the star so Eq. I|19|l is not merely a Poisson type equation 
for It therefore deserves a special treatment (see Ap- 
pendix B in [2(j for a discussion). In the case of strange 
stars, the coefficient = 1/3 in whole star so we have 
to deal with a usual Poisson equation and consequently 
we have to impose a boundary condition for the velocity 
potential at the stellar surface. 

We can define the surface of the star by Jii sur f = uq = 
constant. The surface of the fluid ball is obviously Lie- 
dragged along the fluid 4-velocity vector u, so that this 
last condition gives 



(£ u n)\ s 



= 0, 



(23) 



where £ u is the Lie derivative along the vector field u. 
Let us decompose u in a part along the helical Killing 
vector £ and a part S parrallel to the hypersurface S t : 



u = X(£ + S). 
The condition l|23|l is then equivalent to 
\{£ t n + £ S n)\ sm{ = 0. 



(24) 



(25) 



Now, if the fluid flows obeys to the helical symmetry 
£gn = 0; inserting this relation into Eq. I|25|) leads to 
(£sn)\ BUI { = or equivalently (since S is a spatial vec- 
tor): 



= 0. 



(26) 



Now, let us express S in terms of the spatial vectors U 
and B. First, Eq. 124|) implies n ■ u = An • £. Secondly, 
the fluid motion u can be described by the orthogonal 
decomposition u = T n (n + U) which yields n u = —T n . 
Finally, from Eq. Q, we have n ■ t = —N so that the 
factor A can be expressed as A = T n /N and Eq. I|24() 
becomes 



N 



(£ + S) 



(20) Now, combining Eq. I|27|) and Eq. ©, we have 



u = r„ 



n 



(27) 



(28) 



Comparing with the orthogonal decomposition u = 
T n (n + U), we deduce that S = NU + B. Inserting this 
relation into Eq. (|26|l leads to the boundary condition 



surf 



= 0. 



(29) 



Now, using Eq. (|22|) . we obtain a Neumann- like bound- 
ary condition for $: 



(V nV,*) 



surf 



T n hA 2 



N 



(30) 



surf 
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Considering the elliptic equation (|19fl for "J we see that 
the boundary condition we have obtained is consistent 
with the case n = (or equivalently (H = 0) at the 

stellar surface since, from Eq. H2()|l . V H = ^-V n. 



V. NUMERICAL RESULTS 
A. The method 

The resolution of the above nonlinear elliptic equa- 
tions is performed thanks to a numerical code based on 
multidomain spectral methods and constructed upon the 
Lorene C++ library The detailed description of 

the whole algorithm, as well as numerous tests of the 
code can be found in [46j. Additional tests have been 
presented in Sec. 3 of Q. The code has already been 
used successfully for calculating the final phase of inspi- 
ral of binar y n eutron stars described bypolytropic EOS 
H [H El H! H3 and realistic EOS 0. It is worth 
to stress that the adaptation of the domains (numerical 
grids) to the stellar surface (surface-fitted coordinates) 
used in this code is particulary usefull here, due to the 
strong discontinuity of the density field at the surface of 
strange stars (cf. the top panel in Fig. Adapting 
the grids to the stellar surface allows to avoid the severe 
Gibbs phenomenon that such a discontinuity would nec- 
essary generate when performing polynomial expansions 
of the fields [13. 

The hydrodynamical part of the code has been 
amended for the present purpose, namely to solve 
Eq. (|19|l for the velocity potential ^ subject to the 
boundary condition (|30|) . Let us recall that in the orig- 
inal version of the code, the treatment of Eq. I|19(l was 
different due to the vanishing of the density field at the 
stellar surface (see Appendix B of Ref. |46|h 

We have used one numerical domain for each star and 
3 (resp. 4) domains for the space around them for a small 
(resp. large) separation. In each domain, the number of 
collocation points of the spectral method is chosen to be 
N r x N e x N v = 25 x 17x 16, where N r ,N g , and N v denote 
the number of collocation points (— number of polyno- 
mials used in the spectral method) in the radial, polar, 
and azimuthal directions respectively. The accuracy of 
the computed relativistic models has been estimated us- 
ing a relativistic generalization of the virial theorem j3|| 
(see also Sec. III. A of Ref. Q). The virial relative error 
is a few times 10 -5 for the closest configurations. 



B. Evolutionary sequences 

For each EOS we construct an evolutionary sequence, 
i.e. a sequence of quasi-equilibrium configurations with 
fixed baryon mass and decreasing separation. Such a 
sequence is expected to approximate pretty well the true 
evolution of binary neutron stars, which is entirely driven 
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FIG. 3: Binding energy as a function of gravitational wave 
frequency along evolutionary sequences of corotating bina- 
ries. The solid line denotes strange quark stars, the dashed 
one neutron stars with polytropic EOS, and the dotted one 
point- mass binaries in the 3PN approximation [53^ . The di- 
amonds locate the minimum of the curves, corresponding to 
the innermost stable circular orbit; configurations to the right 
of the diamond are securaly unstable. 



by the reaction to gravitational radiation and hence occur 
at fixed baryon number and fluid circulation. 

For a given rotational state we calculate evolutionary 
sequences of binary system composed of two identical 
neutron stars or two identical strange stars. The evo- 
lution of inspiraling corotating (irrotational) binaries is 
shown in Fig. |3] (Fig. 0J. Fig. and upper panel of 
Fig. 0]show the binding energy -Ebmd versus frequency of 
gravitational waves /gw and lower panel of Fig. 01 show 
the total angular momentum of the systems as a function 
of /gw- The binding energy is defined as the difference 
between the actual ADM mass of the system, Madm, 
and the ADM mass at infinite separation (2.7 Mq in our 
case). The frequency of gravitational waves is twice the 
orbital frequency, since it corresponds to the frequency of 
the dominant part 1 = 2, m = ±2. Solid and dashed lines 
denote quasi-equilibrium sequences of strange quark stars 
binaries and neutron stars binaries respectively. Dotted 
lines in Fig. and Fig. 01 correspond to the 3rd PN 
approximation for point masses derived by [Hflj . Finally 
in Fig. [7|we compare our results with third order post- 
Newtonian results for point-mass particles obtained in 
the effective one body approach by Damour et al. 2000 
[Hij . Damour et al. 2002 [^J and in the standard non- 
resummed post-Newtonian framework by Blanchet 2002 

A turning point of -Ebind along an evolutionary se- 
quence indicates an orbital instability [3(|. This insta- 
bility originates both from relativistic effects (the well- 
known r = 6M last stable orbit of Schwarzschild metric) 
and hydrodynamical effects (for instance, such an insta- 
bility exists for sufficiently stiff EOS in the Newtonian 
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regime, see e.g. |4|| and references therein). It is secular 
for synchronized systems and dynamical for irrotational 
ones. 

In the case where no turning point of -E^ind occurs 
along the sequence, the mass-shedding limit (Roche lobe 
overflow) marks the end of the inspiral phase of the 
binary system, since recent dynamical calculations for 
7 = 2 polytrope have shown that the time to coalescence 
was shorter than one orbital period for configurations at 
the mass-shedding limit [T^J, |2S) ■ Thus the physical in- 
spiral of binary compact stars terminates by either the 
orbital instability (turning point of £"bind) or the mass- 
shedding limit. In both cases, this defines the inner- 
most stable circular orbit (IS CO). The orbital frequency 
at the ISCO is a potentially observable parameter by the 
gravitational wave detectors, and thus a very interesting 
quantity. 
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C. Corotating binaries 

Quasi-equilibrium sequences of equal mass corotating 
binary neutron stars and strange stars are presented in 
Fig. [3J For both sequences we find a minimum of the 
binding energy. In the present rotation state, this lo- 
cates a secular instability |3t| . The important difference 
between neutron stars and strange stars is the frequency 
at which this instability appears. Indeed, there is a dif- 
ference of more than 100 Hz: 1020 Hz for strange stars 
and 1140 Hz for neutron stars. The binding energy is 
the total energy of gravitational waves emitted by the 
system: a corotating binary strange star system emits 
less energy in gravitational waves and loses less angular 
momentum before the ISCO than a binary neutron star 
one with the same mass and compaction parameter in 
infinite separation. 

Comparison of our numerical results with 3rd order 
PN calculations reveals a good agreement for small fre- 
quencies (large separations) (see Fig. |3]and|7|). The 
deviation from PN curves at higher frequencies (smaller 
separation) is due to hydrodynamical effects, which are 
not taken into account in the PN approach. 



D. Irrotational binaries 

In Fig. 0] we present the evolution of the binding en- 
ergy and angular momentum for irrotational sequences 
of binary neutron stars and strange stars. We also verify 
that these sequences are in a good agreement with PN 
calculations for large separations. 

We note important differences in the evolution of bi- 
nary systems consisting of strange stars or neutron stars. 
The strange star sequence shows a minimum of the bind- 
ing energy at /aw — 1390 Hz, which locates a dynamical 
instability 36] and thus defines the ISCO. The minimum 
of i?bind coincides with the minimum of total angular 
momentum J. This is in accordance with the "first law 




1000 1200 



FIG. 4: Binding energy (top panel) and angular momentum 
(bottom panel) as a function of gravitational wave frequency 
along evolutionary sequences of irrotational binaries. The 
solid line denotes strange quark stars, the dashed one poly- 
tropic neutron stars, and the dotted one point-mass binaries 
in the 3PN approximation 53]. The diamonds correspond to 
dynamical orbital instability (the ISCO). 

of binary relativistic star thermodynamics" within the 
IWM approximation as derived by Friedman, Uryu and 
Shibata |3jj and which states that, along an evolutionary 
sequence, 

SM ABM = (31) 

The surface of strange stars at the ISCO is smooth (see 
Fig. EJ. On the contrary the neutron star sequence does 
not present any turning point of -Ebmd 7 so that the ISCO 
in this case corresponds to the mass-shedding limit (final 
point on the dashed curves in Fig. HJ. The gravitational 
wave frequency at the ISCO is much lower for neutron 
star binaries than for strange star binaries. 

As already mentioned the adiabatic index in the outer 
layers of a compact star in a binary system plays a cru- 
cial role in its evolution, especially in setting the mass- 
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FIG. 5: Internal velocity fields of irrotational strange quark 
stars binaries at the ISCO. upper panel: velocity U in the 
orbital plane with respect to the "inertial" frame (Eulerian 
observer); lower panel: velocity field with respect to the coro- 
tating frame. The thick solid lines denote the surface of each 
star. 



shedding limit. Although the crust of a 1.35 M Q neutron 
star contains only a few percent of the stellar mass, this 
region is easily deformed under the action of the tidal 
forces resulting from the gravitational field produced by 
the companion star. The end of inspiral phase of binary 
stars strongly depends on the stiffness of matter in this 
region. It has been shown that the turning-point orbital 
instability for irrotational polytropic neutron stars bina- 
ries can be found only if 7 > 2.5 and if the compaction 
parameter is smaller than certain value (Q, |2{|). In fact, 
as shown in Fig. 31 of paper Q, they didn't find ISCO for 
irrotational binary neutron stars with 7 = 2. 5 or 7 = 3 
for compaction parameter as high as M/R = 0.187. 

In Fig. we present the evolution of two different stel- 
lar radii: the equatorial radius R x , defined as half the 
diameter in the direction of the companion and the po- 
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FIG. 6: Coordinate "radius" (half the coordinate size of a 
star in fixed direction) versus coordinate separation for ir- 
rotational quasi-equilibrium sequences of binary strange stars 
(solid line) and neutron stars (dashed line) . Upper lines corre- 
spond to equatorial radius R x (radius along the x axis going 
through the centers of stars in a binary system) and lower 
lines are polar radius R z (radius along the rotation axis). 



lar radius, defined as half the diameter parallel to the 
rotation axis. For spherical stars R x — R z . We see 
that at the end of the inspiral phase, neutron stars are, 
for the same separation, more oblate (more deformed) 
than strange stars. Binary neutron stars reach the mass- 
shedding limit (the point at which they start to exchange 
matter - a cusp form at the stellar surface in the direction 
of the companion) at coordinate separation d ~ 25km. 
We don't see any cusps for strange stars even for dis- 
tances slightly smaller that the distance corresponding 
to the ISCO - 23.5 km. 

It is worth to remind here the results on rapidly rotat- 
ing strange stars and neutron stars. The Keplerian limit 
is obtained for higher oblatness (more deformed stars), 
measured for example by the ratio of polar and equato- 
rial radius, in the case of stran ge s tars than in the case 
of neutron stars HHHHliiliii^. 

The differences in the evolution of binary (or rotating) 
strange stars and neutron stars stem from the fact that 
strange stars are principally bound by another force than 
gravitation: the strong interaction between quarks. 

As already mentioned the frequency of gravitational 
waves is one of potentially observable parameters by the 
gravitational wave detectors. We can see from Fig. 0that 
the 3rd PN approximations for point masses derived by 
different authors are giving ISCO at very high frequen- 
cies of gravitational waves > 2 kHz. Since in the hy- 
drodynamical phase of inspiral the effect of a finite size 
of the star (e.g. tidal forces) is very important we see 
deviation of our numerical results from point-masses cal- 
culations. The frequency of gravitational waves at the 
ISCO strongly depends on equation of state and the ro- 
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tational state. For irrotational equal mass (of 1.35 M© 
at infinite separation) strange stars binaries described by 
the simple MIT bag model this frequency is ~ 1400Hz 
and for neutron stars binaries described by four different 
realistic EOS it is between 800Hz and 1230Hz HQ. 

VI. SUMMARY AND DISCUSSION 

We have computed evolutionary sequences of irrota- 
tional and corotating binary strange stars by keeping the 
baryon mass constant to a value that corresponds to in- 
dividual gravitational masses of 1.35 Mq at infinite sepa- 
ration. The last orbits of inspiraling binary strange stars 
have been studied in the quasi-equilibrium approxima- 
tion and in the framework of Isenberg- Wilson-Mathews 
approximation of general relativity. In order to calculate 
hydrodynamical phase of inspiraling irrotational strange 
stars binaries, i.e. assuming that the fluid has zero vor- 
ticity in the incrtial frame, we found the boundary condi- 
tion for the velocity potential. This boundary condition 
is valid for both the case of non- vanishing (e.g. self-bound 
matter) and vanishing density at the stellar surface (neu- 
tron star matter). In our calculations strange stars are 
built by strange quark matter described by the simplest 
MIT bag model (assuming massless and non-interacting 
quarks) . 

We have located the end of each quasi-equilibrium se- 
quence (ISCO), which corresponds to some orbital insta- 
bilities (the dynamical instability for irrotational case or 
the secular one for synchronized case) and determined the 
frequency of gravitational waves at this point. This char- 
acteristic frequency yields important information about 
the equation of state of compact stars and is one of the 
potentially observable parameters by the gravitational 
wave detectors. In addition, the obtained configurations 
provide valuable initial conditions for the merger phase. 
We found the frequency of gravitational waves at the 
ISCO to be ~ 1400 Hz for irrotational strange star bi- 
naries and ~ 1000 Hz for synchronized case. The irrota- 
tional case is more realistic since the viscosity of strange 
star matter is far too low to ensure synchronization dur- 
ing the late stage of the inspiral. For irrotational equal 
mass (of 1.35 A/q) neutron star binaries described by re- 
alistic EOS the frequency of gravitational waves 
at the ISCO is between 800Hz and 1230Hz, much lower 
than for a binary strange quark star built of self-bound 
strange quark matter. We have considered only strange 
quark stars described by the simple MIT bag model with 
massless and non-interacting quarks. In order to be able 
to interpret future gravitational-wave observations cor- 
rectly it is necessary to perform calculations for different 
strange star EOS parameters (taking also into account 
the existence of a thin crust) and for large sample of 
neutron stars described by realistic equations of state. 



For some MIT bag model parameters one is able to ob- 
tain less compact stars than considered in the present 
paper. In this case the frequency of gravitational waves 
at the end of inspiral phase will be lower than obtained 
by us. It should be also taken into account that stars in 
a binary system can have different masses j^Oj • The case 
of binary stars (with equal masses and different masses) 
described by different strange quark matter models will 
be presented in a separate paper [6l| . 

We have shown the differences in the inspiral phase 
between strange quark stars and neutron stars described 
by polytropic equation of state having the same gravita- 
tional mass and radius in the infinite separation. It was 
already shown by Bejger et al. 2005 that the frequency 
of gravitational waves at the end point of inspiraling neu- 
tron stars described by several realistic EOS without ex- 
otic phases (such as meson condensates or quark matter) 
can be predicted, in a good approximation, by studying 
binaries with assumed polytropic EOSs with 7 = 2 or 2.5. 
For realistic EOS and polytropes with 7 < 2.5 pi |25| 
a quasi-equilibrium irrotational sequence terminates by 
mass-shedding limit (where a cusp on the stellar surface 
develops). 

We found that it wasn't the case for inspiraling strange 
star binaries which are self-bound objects having very 
large adiabatic index in the outer layer. For both syn- 
chronized and irrotational configurations, we could al- 
ways find a turning point of binding energy along an 
evolutionary sequence of strange quark stars, which de- 
fines an orbital instability and thus marks the ISCO in 
this case. In the irrotational case for the same separation 
strange stars are less deformed than polytropic neutron 
stars and for the same ratio of coordinate radius R x /R z 
their surfaces are more smooth. A cusp doesn't appear 
on the surface of a strange star in a binary system even 
for separation corresponding to orbital instability. The 
frequency of gravitational waves at the end of inspiral 
phase is higher by 300 Hz for the strange star binary sys- 
tem than for the polytropic neutron star binaries. The 
differences in the evolution of binary (or rotating) strange 
stars and neutron stars stem from the fact that strange 
stars are principally bound by an additional force, strong 
interaction between quarks. 
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FIG. 7: Binding energy versus frequency of gravitational waves along evolutionary sequences of corotational (thick dashed 
lines) and irrotational (thick solid lines) equal mass (of 1.35 Mq) strange stars and polytropic neutron stars binaries compared 
with analytical results at the 3rd post-Newtonian order for point-masses by Damour et al. 2000 |H^|. Damour et al. 2002 |55| 
and Blanchet 2002 [H 
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